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In order to clarify the mechanism of the enhancement of superconducting transition tempera- 
ture T c due to anharmonic local oscillation of a guest ion in a cage composed of host atoms, i.e., 
rattling, we analyze the anharmonic Holstein model by applying the Migdal-Eliashberg theory. 
From the evaluation of the normal-state electron- phonon coupling constant, it is found that the 
strong coupling state is developed, when the bottom of a potential for the guest ion becomes 
wide and flat. Then, T c is enhanced with the increase of the anharmonicity in the potential, 
although T c is rather decreased when the potential becomes a double-well type due to very 
strong anharmonicity. From these results, we propose a scenario of anharmonicity-controlled 
strong-coupling tendency for superconductivity induced by rattling. We briefly discuss possible 
relevance of the present scenario with superconductivity in /3-pyrochlore oxides. 
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1. Introduction 

Recently, materials with cage structure have attracted 
much attention in the research field of condensed mat- 
ter physics from the emergence of exotic magnetism 
and novel superconductivity. The materials are, for in- 
stance, filled skutteruditcs, 1-3 ' /3-pyrochlore oxides, 4-7 ) 
and clathrate compounds. 8 ~ 10 ' A common feature of this 
material group is the existence of nano-size cage com- 
posed of relatively light atoms. The ion contained in the 
cage, frequently called a guest ion, feels a highly anhar- 
monic potential and it can oscillate with large amplitude. 
Such oscillation of the guest ion in the cage is called 
rattling, which is expected to be an origin of interesting 
physical properties of cage materials such as magnetically 
robust heavy electron state of filled skutterudites, 11 ' su- 
perconductivity of /3-pyrochlore oxides, 12 ' and high ther- 
moelectric performance of clathrate compounds. 13 -' 

From a theoretical viewpoint, heavy electron state 
with phonon origin has been actively investigated in 
recent years. The magnetically robust heavy electron 
phenomenon has been discussed from the viewpoint of 
the Kondo effect with phonon origin. 14 ~ 16 ' The periodic 
Anderson- Holstein model has been analyzed with the use 
of a dynamical mean-field theory and a mechanism of the 
mass enhancement due to large lattice fluctuations and 
phonon softening towards double- well potential has been 
addressed. 17 ^ One of the present authors has performed 
numerical calculations for the Anderson model coupled 
with local Jahn- Teller and Holstein phonons. 18-21 ' From 
the evaluation of electronic specific heat coefficient j c in 
the Andcrson-Holstcin model, it has been found that 7 e 
is enhanced by the rattling, namely, the oscillation of 
the guest ion in a potential with wide and flat bottom. 
Then, it has been shown that 7 enhanced by the rattling 
is actually robust for an applied magnetic field. 22 ) Fur- 
thermore, it has been pointed out that the Kondo effect 
due to the rattling exhibits peculiar isotope effect and 
it can be the experimental evidence for rattling-induced 



heavy fermion phenomenon. 23 ' 

Now we turn our attention to rattling-induced phe- 
nomena in /3-pyrochlore oxides AOS2O6 (A=K, Rb, and 
Cs). It has been observed that T c increases with the de- 
crease of radius of A ion: T C =9.6K for A=K, T C =6.4K 
for A=Rb, and T C =3.25K for A=Cs. 12 ' The origin of this 
difference in T c has been expected to be the anharmonic 
oscillation of A ion. In fact, it has been revealed that 
the anharmonicity of the potential for A ion is enhanced 
when we change A ion in the order of Cs, Rb, and K. 24 ' In 
particular, KOS2O6 has relatively high T c in comparison 
with other /3-pyrochlore oxides and it exhibits a remark- 
able first order transition without symmetry change at 
T p = 7.5if . 25 ' In /3-pyrochlore oxides, characteristic be- 
havior has been observed in the temperature dependence 
of electrical resistivity 26 ' and NMR relaxation rate. 27 ' 

These interesting phenomena in /3-pyrochlore oxides 
have been considered to originate from the rattling, i.e., 
anharmonic oscillation of alkali ion in the cage. Dahm 
and Ueda have explained unusual temperature depen- 
dence of electrical resistivity and NMR relaxation rate 
by anharmonic phonons. 28 ' Hattori and Tsunetsugu have 
investigated more realistic model including three dimen- 
sional anharmonic phonons in tetrahedral symmetry and 
have found that T c is strongly enhanced with increas- 
ing the third-order anharmonicity of the potential. 29,30 ' 
Chang et al. have discussed the superconductivity and 
the first-order transition in KOS2O6 by using the strong- 
coupling Eliashberg approach in the anharmonic Hamil- 
tonian including forth-order terms. 31 ' As for the first- 
order transition observed in KOS2O6, Fuse and Ono have 
proposed a possible scenario on the basis of bipolaron for- 
mation. 32 ' From these efforts, it has been gradually rec- 
ognized that the anharmonic oscillation of the ion in the 
cage is closely related to the emergence of superconduc- 
tivity with relatively high T c . However, the mechanism 
of the enhancement of T c due to anharmonicity is still 
unclear. In particular, it is worth while investigating the 
relation among T c , mass enhancement, and anharmonic- 
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ity in the same model Hamiltonian. 

In this paper, we analyze the anharmonic Holstein 
model in an adiabatic approximation, in order to un- 
derstand the enhancement of superconducting transition 
temperature T c and the electron mass enhancement. For 
the purpose, we calculate the electron self-energy by us- 
ing standard Green's function method and evaluate T c 
by solving the Eliashbcrg equation within the Migdal ap- 
proximation. 33,34 ) Due to the evaluation of an effective 
electron-phonon coupling constant A n , it is found that A n 
is increased with the increase of the anharmonicity even 
if the original electron-phonon coupling is weak. This 
strong-coupling tendency due to anharmonicity induces 
the enhancement of electron effective mass and T c . Then, 
it is shown that T c becomes maximum when the guest ion 
oscillates with large amplitude in a potential with a wide 
and flat bottom. We also find that the magnitude of T c 
is increased with the increase of the Debye- Waller factor, 
as observed experimentally in /3-pyrochlore oxides. 

The organization of this paper is as follows. In Sec. 2, 
we introduce the model Hamiltonian and discuss briefly 
the relation between potential shape and anharmonicity. 
Then, we provide the formulation to evaluate electron 
self-energy and T c within a framework of the Migdal- 
Eliashberg theory. In Sec. 3, we show our calculated re- 
sults and discuss the key role of anharmonicity for the 
emergence of rattling-induced superconductivity. We also 
discuss possible relevance of our scenario with experi- 
mental results of /3-pyrochlore oxides. Finally, we sum- 
marize this paper in Sec. 4. Throughout this paper, we 
use such units as h=kB=l. 

2. Model and Formulation 

2.1 Anharmonic Holstein Model 

Let us introduce the anharmonic Holstein model in 
which conduction electrons are coupled with local os- 
cillation of guest ion in the anharmonic potential. The 
Hamiltonian consists of three parts as 

H = £kc{ a c ka + ff eph + iJp h , (1) 

where e k is the energy of conduction electron, k is mo- 
mentum, Cfco- indicates an annihilation operator of elec- 
tron with spin a and momentum k, iJ cp h indicates the 
electron-phonon coupling, and H p h is the phonon part. 
Throughout this paper, we consider a half-filling situa- 
tion by appropriately adjusting the value of a chemical 
potential included in although we do not explicitly 
include the chemical potential term in H . 

The second term -ff op h in eq. (1) indicates the electron- 
phonon coupling part, given by 



H, 



eph 



(2) 



where g denotes an electron-phonon coupling constant, 
i indicates an ion site, c\„ is the annihilation operator 
of electron at a site i, and Q\ is normal coordinate of 
oscillation of the guest ion at a site i. Note here that the 
reduced mass for the normal oscillation is set as unity in 
this paper for simplicity. 



The third term H v \ x in cq. (1) denotes the phonon part, 
written as 



H, 



ph 



[lf/2 + ViQi)], 



(3) 



where Pi denotes the canonical momentum of the guest 
ion at a site i and V(Qi) indicates the potential for the 
guest ion. In order to include the effect of anharmonicity, 
we express V(Qi) as 

V(Qi) = LolQl/2 + k A Qt + keQt (4) 

where luq is the energy of the guest ion and k^ and k§ 
denote the coefficients for fourth- and sixth-order anhar- 
monic terms, respectively. 

For quantum statistical calculations, it is convenient 
to introduce the annihilation operator of phonon through 
the standard procedure as Qi = (oti +a\ )/ yj (2wo). Then, 
we obtain i/ cp h and £f p h, respectively, as 



(5) 



and 



H ph =u J0 J2i a l a i + 1 /' 2 +f 3 ( a i + «i) 4 +7(ai + 4)% ( 6 ) 



where a, j3, and 7 are given by 



8< 



respectively. Note that a is non-dimensional electron- 
phonon coupling constant, while (3 and 7 indicate non- 
dimensional fourth- and sixth-order anharmonicity pa- 
rameters, respectively. By using these non-dimensional 
parameters, we express the potential as 



V{qi) = auj {q? + IQapqf + 64a 2 7 gf ), 



(8) 



where qi = Qi/£ and t denotes a typical length scale of 
the oscillation, given by I — g/uiQ. 

Before proceeding to the formulation, here we briefly 
discuss how the potential shape is changed by anhar- 
monicity parameters. It is noted that we consider only 
the case of 7 > 0, since the oscillation of the guest ion 
should be confined in a finite space of the cage. In or- 
der to characterize the potential shape, we define the 
re-scaled anharmonicity parameter as 



(9) 



With the use of this parameter, even if the value of 
7 is changed, the potential shape can be discussed in 
the same range of f3' . As already mentioned in previous 
papers, 22 ' 23 ' 35 -' the potential shapes are classified into 
three types by the values of as on-center type for 
P' > -v / 3/2, rattling type for -y/3/2 > /3' > -1, and 
off-center type for /3' < — 1. Typical potential shapes are 
shown in Fig. 1 for various values of /?'. Note that the 
potential height Vq is defined as the difference between 
the potential values at local maxima and minima. As a 
typical example, we show Vb for the case of j3' = —y/3/2. 
Here we note that (i' = — y3/2 is defined by the value 
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consider an adiabatic situation in which W is much larger 
than the phonon energy l)q. 

Concerning the phonon Green's function Dq of anhar- 
monic system, it is convenient to express it in the spectral 
representation as 



Pph(w) 



(12) 



where v n is the boson Matsubara frequency defined by 
v n = 2mrT, and p p h(u)) is phonon spectral function, 



given by 



Pph(w) 



Aij6(uj + Ei- Ej). 



(13) 



-4 



-2 







Here Ei is the i-th eigenenergy of £f p h and the spectral 
weight Aij is given by 



1 



Fig. 1. (Color online) Anharmonic potentials for 7 = 10~ 5 . We 
show the harmonic potential as a dotted curve for comparison. 
The curve of /3' = —0.01 denotes the on-center type, while that 
of 0' = —1.2 indicates the off-center one. As for the meaning of 
Vb, see the main text. 



at which the potential exhibits saddle points at (ft w ±2, 
while /3' = — 1 indicates the case for which all the poten- 
tial values at three minima become zeros. It is observed 
that the rattling-type potential has relatively wide and 
flat region in the bottom in comparison with other two 
cases. Namely, in this case, the guest ion is considered to 
oscillate with large amplitude in the anharmonic poten- 
tial. From the purpose of this paper, it is quite natural 
to focus on the parameters around the values for the 
rattling-type potential. 

2.2 Formulation 

In order to discuss the superconducting transition tem- 
perature T c , we analyze the Hamiltonian by applying the 
formalism of the Migdal-Eliashbcrg theory 33, 34 ) First we 
consider the normal self-energy in the Migdal approxi- 
mation in which electron-phonon vertex corrections are 
ignored in an adiabatic limit of loq -C W, where ojq de- 
notes the phonon energy and W is the electron band- 
width. Then, in the second-order perturbation theory in 
terms of g, the self-energy £ is given by 



(14) 



where is the i-th eigenstate of H p h, il is the partition 
function, given by Q — y^ - e~ Ei / T , and m is the non- 
dimensional displacement operator, given by ui=ai + a\. 

Let us define an electron mass enhancement factor 
given by 

Y,{inT) 



Z = 1 - 



inT 



(15) 



We note that the effective mass m* is given by m*=Zm, 
where m is bare electron mass. Here we assume that Z 
is given by the difference of £ on the imaginary axis at 
the lowest temperature, although Z should be defined by 
the energy differentiation of E on the real axis as 

<9Re£H 



Z= 1- 



(16) 



j=0 



£(«w„) 



where S(w) denotes the electron self-energy on the real 
axis and w denotes the energy. Since S(w) is found to 
exhibit a peak in the order of the phonon energy ujq, we 
replace the energy differentiation with the energy differ- 
ence of the width T for T < wo- 35 ^ In this paper, we are 
interested in the quantity around T = T c and T c is con- 
sidered to be lower than luq. Thus, we adopt the present 
definition for Z. 

For the purpose to evaluate T c , we solve the linearized 
gap equation at T=T C , given by 

^EE D oK-iw„:)G(fc',^), (10) <p{iu n )=aulTY^^D {iuj n -iu n ,)F{k',iu n ,), (17) 



where T is a temperature, uj n is the fermion Matsubara 
frequency defined by u> n = (2n + l)wT with an integer 
n, Dq is the phonon Green's function shown below, and 
G(fc, iu) n ) is the electron Green's function, given by 

G{k, iu„) = -. 1 w . . ■ (11) 

lUJ n - £fe - M ?w nJ 

We note that the site dependence dose not appear in 
Dq, since we consider Einstein- type local phonons in the 
adiabatic approximation. Thus, the self-energy does not 
depend on the momentum. Note also that in our cal- 
culation, we use the bare phonon Green's function Do, 
instead of dressed phonon Green's function D, since we 



where <fi(iuj n ) is anomalous self-energy and F(k,iuj n ) is 
anomalous Green's function, given in the linearized form 
in the vicinity of T c as 



F(k,iuj n ) = -G{k,iuj n )G{-k,-iuj n )4)(iuj n ). 



(18) 



Note that here we assume the isotropic s-wave gap, since 
we consider Cooper pairs mediated by phonons. Thus, 
the momentum dependence is suppressed in the anoma- 
lous self-energy. 

The calculation procedure is as follows: First we deter- 
mine the normal self-energy by solving sclf-consistently 
cqs. (10) and (11). Then, by using the obtained normal 
self-energy, we solve the gap equation eqs. (17) and (18). 
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Note that T c is determined as a temperature at which 
the positive maximum eigenvalue of the gap equation 
becomes unity. 

In the actual calculations, we assume the electron den- 
sity of states with rectangular shape of the electron band- 
width W. Hereafter W is taken as the energy unit and 
we set W = 1. For the sum on the imaginary axis, we 
use 32768 Matsubara frequencies. In order to accelerate 
the actual calculations in eqs. (10) and (17), we exploit 
a Fast-Fourier- Transformation algorithm for the summa- 
tion in terms of n' . The eigenvalue of the gap equation 
is obtained with the use of the power method. For diag- 
onalization of -ff pn , we use 250 phonon basis. Note that 
we have checked the convergence of the eigenenergy by 
comparing the results obtained by using 300 phonon ba- 
sis. Throughout this paper, we set luq = 0.1 and A = 0.5, 
where A is the Eliashberg electron-phonon coupling con- 
stant given by A = 2auo/W. 

3. Calculated Results 

3.1 Superconducting transition temperature 

Let us first discuss the effect of anharmonicity on the 
superconducting transition temperature T c . In Fig. 2(a), 
we depict three curves of T c vs. ft for 7=0.01, 0.001, 
and 0.0001. For 7=0.0001, T c slowly increases with the 
decrease of ft in the range of ft > —1 and it rapidly de- 
creases at ft = — 1. Note that we could not calculate T c 
smaller than 0.001 due to the limitation of the present 
numerical calculations. For 7 = 0.001, we also observe 
that T c is enhanced by the increase of fourth-order an- 
harmonicity and it takes the largest value at ft ~ — 1. 
Note that the largest value is also the highest among 
those in three curves. In the parameter region of the off- 
center type potential, T c rapidly decreases. 

For 7 = 0.01, T c also exhibits a peak structure, but 
the width of the peak is broad and T c becomes highest 
for ft « —1.2. As understood from Fig. 1, the potential 
is the off-center type with a couple of shallow minima 
at q\ f=a ±2. Although T c is not always maximum just 
in the range of -1 < ft < -y/3/2, we conclude that 
the highest T c is obtained in the vicinity of the rattling 
potential with a wide and flat bottom. 

In order to consider further the mechanism of the in- 
crease of T c , we estimate the effective electron-phonon 
coupling constant A n at T = T c . With the use of the 
phonon spectral function p p h(w), A n is evaluated as 

PphH 



_ \uj 



-dui. 



(19) 



For the case of /3=7=0 (harmonic phonon), we easily 
obtain A n =A from the definition of A n and the phonon 
spectral function. Note, however, that in general, A n is 
different from A. Thus, the difference between A n and A 
is considered to indicate the effect of anharmonicity on 
the electron-phonon coupling constant. 

The calculated results for A n are shown in Fig. 2(b). 
We immediately observe that in the on-center type po- 
tential {ft > —a/3/2), A n is not so enhanced from the 
bare value of A. When ft enters the value of the region 
of the rattling-type potential (— a/3/2 > ft > —1), A n is 
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Fig. 2. (Color online) (a) Superconducting transition tempera- 
ture T c vs. p' for 7=0.01, 0.001, and 0.0001. Here we set A = 0.5 
and wo = 0.1. (b) Effective electron-phonon coupling constant 
A n vs. f) 1 for the same parameters as in the panel (a). 



rapidly increased with the decrease of ft. When we fur- 
ther decrease the values of ft, A n continues to increase 
even in the off-center type potential (ft < —1). Here we 
should pay due attention to the validity of the present 
adiabatic approximation in the off-center type potential, 
but this point will be commented later. At the present 
stage, we suggest that the enhancement of T c is related 
to the strong coupling state due to anharmonicity. 

Concerning the formation of the peak in T c , we provide 
a qualitative comment from the strong-coupling effect. In 
the BCS theory, it is well known that T c is expressed as 
T c = 1.13cj e- 1/A . When we consider the renormalization 
effect, the formula of T c is changed as T c = 1.13u)oe~ z ^ n ■ 
where Z is the mass enhancement factor. Namely, l/Z 
indicates the renormalization effect. Note that when A n is 
moderately large, Z is roughly given by 1 + A n , leading to 
the simplest version of the McMillan formula. 36 ' When 
A n is small, the renormalization effect is not so significant 
and the increase of A n contributes to the enhancement 
of T c . On the other hand, when A n becomes much larger 
than unity, the renormalization effect should be domi- 
nant and it rather suppresses the value of T c . Thus, the 
peak is formed in T c as a function of ft. 

3.2 Phonon softening due to anharmonicity 

In the previous subsection, it has been clarified that 
the enhancement of T c is triggered by the strong-coupling 
tendency due to anharmonicity. Next we attempt to un- 
derstand the origin of such tendency For the purpose, 
let us examine the change of the phonon spectral func- 
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Fig. 3. (Color online) Phonon spectral function p p h(w) vs. energy 
w for (a) 0' = -0.1, (b) 0' = 0.88, (c) 0' = -1.0, and (d) 
0' = —1.1. Other parameters are set as A = 0.5, ujq = 0.1, and 
7 = 0.001. Note that we calculate p p ^(uj) at T = T c . 



tion, when we change the value of /?'. In Figs. 3, we 
depict Pph(w) for /3'=-0.1, -0.88, -1.0, and -1.1 at 
T = T c with A=0.5, w=0.1, and 7=0.001. The panel (a) 
for f3' = —0.1 indicates the results for the on-center type 
potential. We observe that the phonon spectral function 
is quite similar to that of the harmonic phonons with a 
finite weight at lo=loq. Then, A n for this case is almost 
equal to A, as observed in Fig. 2(b). 

When we decrease the value of /3', the phonon state 
with finite weights are gradually shifted to zero, as ob- 
served in Figs. 3(b) and 3(c), which are the results for the 
rattling-type potential. It is easy to understand that A n 
becomes large, since low-energy phonon states with finite 
weights contribute to A n through p p h(w)/w in eq. (19). 
In Fig. 3(d) for the case of (3' = —1.1, the largest weight 
of pph(w) shifts to the lower energy and the moderate 
weights appear in higher u>. Thus, A n becomes large in 
the parameter region of the off-center type potential. In 
short, the increase of the anharmonicity leads to phonon 
softening. Namely, the weight of p p h(w) is shifted to lower 
energy and thus, A n is enhanced in the anharmonic po- 
tential. Note that our conclusion is consistent with the 
previous result. 17 ^ 

Next we discuss the phonon excitation in order to un- 
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Fig. 4. (Color online) Five lowest eigenenergies Eo ~ E4 (solid 
curves with symbols) and Vq vs. 0' (broken curve) for (a) 7=0.01, 
(b) 0.001, and (c) 0.0001. Note that we do not show all eigenener- 
gies in these figures. We set A = 0.5 and ujq = 0.1. Since Vo dose 
not exist for 0' > — v3/2 due to its definition, Vo is depicted in 
the region of 0' < -yU/2. 



derstand the relation between T c and anharmonic oscilla- 
tion. In Figs. 4(a)-4(c), we depict eigenenergies Eq ~ E4 
ofH ph and Vo for 7 =0.0001, 0.001, and 0.01, respectively 
First let us consider overall tendency of the change of 
eigenenergy due to the anharmonicity found in Figs. 4(a)- 
4(c). In the case of the on-center type potential for 
fi' > — V3/2, all eigenenergies are gradually decreased 
with the decrease of /?'. For -1.0 < 0' < -\/3/2, in 
which we observe the rattling-type potential, the width 
between eigenenergy levels becomes small. This is the ori- 
gin of the shift of phonon spectral function, as observed 
in Fig. 3. For /?' < —1.0, the eigenenergies become doubly 
degenerate, which is originating from the left and right 
deep minima in the off-center type potential. The width 
between degenerate energy levels are rather large, when 
we further decrease /?', but the degenerate phonon states 
lead to the strong coupling state in the present adiabatic 
approximation. This point will be commented later. 

Intuitively, the phonon softening due to anharmonic- 
ity is understood from the energy of oscillator confined in 
a potential. If the potential size becomes wide, the mo- 
mentum of the oscillator becomes small and the kinetic 
energy is suppressed due to the well-known uncertainty 
principle. In the present case, the size of the potential 
bottom is effectively enlarged due to the anharmonicity 
and thus, the typical phonon energy is suppressed due to 
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the quantum effect. We believe that it is the main effect 
of anharmonicity on the electron-phonon coupling state. 

Next we turn our attention to the change of Vo, since 
Vq competes with the zero-point energy Eq. Note again 
that Vo is defined as the difference between the potential 
values at local maxima and minima. From the compe- 
tition between Vq and Eq, we can discuss the quantum 
effect of potential shape, when there exists structure in 
the bottom of the potential. For 7=0.0001 in Fig. 4(a), 
since Vq is always larger than Eq and E\, the guest ion 
oscillates in the narrow width of the potential bottom, 
which is almost the same as that of the harmonic poten- 
tial in the low-temperature region. On the other hand, 
for /?' < —1.0, since the phonon energies are degener- 
ate, the oscillation state is drastically changed from the 
small oscillation at the potential minimum to the quan- 
tum tunneling or thermal activation between two poten- 
tial minima. Thus, the behavior of T c is rapidly changed 
at j3'=— 1.0 in the present calculations. For 7=0.001 in 
Fig. 4(b), Vq is close to E at /?' ~ —1.0, corresponding 
to the peak observed in the curve for T c in Fig. 2(a). For 
7=0.01, Vq crosses E at /3' ~ —1-2, which also corre- 
sponds to the peak position in T c . 

From these results, we consider that the behavior of 
T c is closely related to E and Vq. For 7=0.001 and 0.01, 
it is found that T c takes the maximum value at the pa- 
rameters corresponding to So ~ Vq. In this case, the 
oscillation amplitude of guest ion is rapidly increased, 
since the ion can feel the potential structure in the en- 
ergy region less than the zero-point energy. Namely, the 
condition of Eq ~ Vo just indicates the rattling oscilla- 
tion of the guest ion. Then, T c is considered to take the 
maximum value in the rattling type potential. 

Hattori and Tsunetsugu have discussed the variation 
of T c by using three-dimensional anharmonic phonons 
in the tetrahedral symmetry. 30 -' They have found a 
crossover in the energy spectrum to the quantum tunnel- 
ing regime and T c becomes large at the crossover point. 
We consider that their condition of the crossover point 
corresponds to Eq ~ Vo- Thus, our result is consistent 
with the previous study. 

3.3 Zero-point energy and anharmonicity 

Now we focus on the relation between the zero-point 
energy Eq and the potential height Vq- When the local 
maximum dose not exist for fi' = — \/3/2, Vo is defined by 
the potential height at the saddle point. For —y/3/2 < (3' , 
Vo cannot be defined. It should be noted that Eq is the 
lowest energy for the oscillation of the guest ion. 

For the rattling-type potential, the potential has min- 
imum at qi = 0. At low temperatures, for Eq < Vo, the 
guest ion oscillates in the narrow width of the potential 
bottom, which is almost the same for that of the har- 
monic potential. On the other hand, for E > Vo, the 
guest ion can oscillate with large amplitude, since the 
potential is wide for the energy larger than Vo- Also in 
the off-center type potential, for Eq > Vo, the guest ion 
can oscillate with large amplitude, since the structure at 
the bottom is masked by the zero-point energy. However, 
for Eq < Vq, the guest ion oscillates in one of potential 
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Fig. 5. (Color online) (a) Superconducting transition tempera- 
ture T c vs. 7. (b) Mass enhancement factor Z vs. 7. (c) Effec- 
tive electron-phonon coupling constant A n at T = T c vs. 7. (d) 
Non-dimensional Dcbyc- Waller factor / vs. 7. (e) The zero-point 
energy Eq and the potential height at the saddle point Vo vs. 7. 
For these panels, we always keep the relation of /?' = — >/3/2. 
Other parameters are set as A=0.5 and wo=0.1. 

minima and sometimes moves to other potential minima 
due to quantum tunneling effect. Thus, it is considered 
that the oscillation of the guest ion in the anharmonic 
potential is quite different between two cases of E > Vo 
and E <V . 

Here let us discuss T c by focusing on the potential at 
P' = — v / 3/2, which is depicted by thick curve in Fig. 1. 
Note again that Vo is the potential height at the sad- 
dle point in this case and (3 depends on 7 through the 
relation of (3 = — ^/3j/2. In Figs. 5(a)-5(c), we depict 
7 dependence of T c , Z, and A n , respectively. When we 
increase the value of 7 from 10 -5 , all those values are 
increased, but they turn to be decreased to form peaks 
at a certain value of 7, which is defined as 7*. 
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In order to elucidate the relation with the amplitude 
of the guest ion, we evaluate the Dcbye- Waller factor 
which represents the intensity of thermal motion. In the 
non-dimensional form, it is given by 

_FW = M1 (20 ) 
J \G\ 2 £ 2 3£ 2 ' v ' 

where Fdw denotes the Debye- Waller factor, G denotes 
the reciprocal lattice vector, £ is the typical length scale 
defined by £=g/u> 2 , and (■ • • ) denotes the operation to 
take thermal average. In Fig. 5(d), we show the result 
for /. We find that / exhibits the peak structure which 
is quite similar to those in Figs. 5(a)-5(c). Namely, / is 
closely related to the behavior of T c , Z, and A n . 

In Fig. 5(e), we depict two curves of Eq and Vq. At 
7 = 7*, they cross with each other, so that the guest 
ion oscillates with large amplitude in the case of Eq ~ 
Vq. Thus, the large amplitude of the guest ion leads to 
the strong coupling state and the enhancement of T c . 
In addition, it is found that this oscillation causes the 
heavy electron state. This is consistent with our previous 
result. 35 ^ 

3-4 D ebye- Waller factor and T c 

Let us now try to understand the tendency of the vari- 
ation of T c experimentally observed in /3-pyrochlore ox- 
ides AOS2O6 (A = K, Rb, Cs) by focusing on the Debye- 
Waller factor at the room temperature and the anhar- 
monicity. 

First we briefly summarize the experimental facts 
for /3-pyrochlore oxides. The isotropic atomic displace- 
ment parameter U- lso at the room temperature which 
corresponds to the Debye- Waller factor in our calcula- 
tion increases as 2.5 x 10 _4 nm 2 , 3.41 x 10~ 4 nm 2 , and 
7.4 x 10 _4 nm 2 in the order of decreasing ionic radius, 
i.e., Cs, Rb, and K. 12 ) On the other hand, T c and A n 
for /3-pyrochlore oxides are summarized as T c =3.25 K 
and A n =2.76 for Cs, T c =6.4 K and A n =3.38 for Rb, and 
T c =9.6 K and A n =6.3 for K. Note that the experimen- 
tal value of A n is obtained by 1 + A n =7cxp/7b, where 
7exp is the experimental Sommerfeld constant and 7b is 
evaluated by the band-structure calculations. Moreover, 
Kuncs et al. have clarified that the anharmonicity of the 
potential which A ion feels increases in the order of Cs, 
Rb, and K. 24 ) Thus, in /3-pyrochlore oxides, T c and A n 
are closely related to the anharmonicity. 

In Figs. 6(a) and 6(b), we show the calculated results of 
T c vs. / and A n vs. /, respectively. Here we change the an- 
harmonicity parameters (/?, 7) on the straight line which 
connects the points of {/3, 7 )=(-0.00392, 0.00154) and 
(—0.00253, 0.00451). Namely, since the non-dimensional 
anharmonicity parameter LV is changed from —0.35 to 
— 1.0, the potential shape is changed from on-center type 
to rattling one. It is understood that when / is increased, 
both T c and A n are increased. This tendency agrees well 
with the experimental situation. Note that it is easy to 
find other parameter regions to depict the figure similar 
to Fig. 6. It is concluded that in the rattling potential 
with large amplitude of the guest ion, T c and the mass 
enhancement is increased with the increase of the Debye- 
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Fig. 6. (Color online) (a) Superconducting transition tempera- 
ture T c vs. / and (b) A n vs. /. Concerning the parameters to 
depict these curves, see the main text. 

Waller factors, as observed in /3-pyrochlore oxides. 

4. Discussion and Summary 

In this paper, we have analyzed the anharmonic Hol- 
stein model which describes the coupling of conduction 
electrons with local ion oscillations in the potential with 
fourth- and sixth-order anharmonic terms, in order to un- 
derstand the enhancement of T c due to anharmonicity. 
First we have investigated T c and the effective electron- 
phonon coupling constant A n . Then, we have found that 
T c is enhanced by the strong coupling tendency due to 
anharmonicity. When the guest ion oscillates with large 
amplitude, the phonon softening occurs. The phonon 
spectral weights are shifted towards lower energy and the 
strong coupling situation is realized even if we consider 
originally the weak-coupling case. This strong-coupling 
tendency due to anharmonicity is considered to be a 
route to arrive at high T c . We believe that this type 
of strong-coupling superconductivity is realized in /?- 
pyrochlore oxides. The relation between T c vs. / as well 
as A n vs. / seems to support our present proposal. 

In the present calculation, wc have always used the 
same phonon energy ujq, but in actuality, u>q is differ- 
ent from material to material. The energy of A ion has 
been reported in the specific heat measurements, 37 ' the 
Raman scattering experiment 38 -* and the neutron scatter- 
ing experiment. 39,40 ) Nagao et al. have reported that the 
energies of the oscillator are 75. IK for A=Cs, 66. 4K for 
A=Rb, and 22K and 61K for A=K in /3-pyrochlore oxides 
AOS2O6. 37 ' If the oscillation energy for the case of A=K 
is taken as 61K, roughly speaking, it is validated to use 
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the same ujq for three /3-pyrochlore oxides. However, if we 
pay our serious attention to the phonon energy of 22K, 
we should consider a situation in which the oscillators 
for A=Cs and Rb have same energy, while the energy for 
A=K is smaller than it. If we attempt to pursue quan- 
titative agreement with the experimental facts, it may 
be necessary to consider the difference among materials. 
However, the research along such a direction should be 
done with the use of more realistic model. It is one of 
future problems. 

In the present adiabatic calculations, A n is monoton- 
ically increased with the decrease of j3', as shown in 
Fig. 2(b). We consider that the calculations are valid for 
the on-center type, the rattling-type, and the off-center 
type potential with shallow minima. However, when the 
value of A n is enhanced in the off-center type poten- 
tial with deep minima at non-zero (fr's, another effect 
should be significant. Namely, for large A n , the effective 
bandwidth given by W/Z becomes comparable with the 
phonon energy and the adiabatic condition will be vi- 
olated. Thus, in the strong-coupling region, it is essen- 
tial to include the non-adiabatic effect through electron- 
phonon vertex corrections. When we include such non- 
adiabatic effect, the degeneracy in phonon energy origi- 
nating from the right and left positions of the off-center 
type potential will be lifted by the electron-phonon cou- 
pling. In the present calculation, unfortunately, the ver- 
tex corrections are not included and the results in the 
region of \/3'\ 1 with negative ft' should be reconsid- 
ered with the inclusion of the vertex corrections. This is 
another future task. 

In summary, we have found that the large ampli- 
tude anharmonic oscillation induces the enhancement 
of T c and the heavy electron state through the strong- 
coupling tendency controlled by the anharmonicity. We 
have shown that when the Debye- Waller factor is in- 
creased, both T c and A n are increased. This tendency 
agrees qualitatively well with the experimental facts in 
/3-pyrochlore oxides. 
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